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Abstract. The onset of thermal convection in a horizontal layer of a colloidal suspension is investigated in 
terms of a continuum model for binary-fluid mixtures where the viscosity depends on the local concentration 
of colloidal particles. With an increasing difference between the viscosity at the warmer and the colder 
boundary the threshold of convection is reduced in the range of positive values of the separation ratio ip 
with the onset of stationary convection as well as in the range of negative values of with an oscillatory 
Hopf bifurcation. Additionally the convection rolls are shifted downwards with respect to the center of the 
horizontal layer for stationary convection {tp > 0) and upwards for the Hopf bifurcation {tp < 0). 

PACS. 47. 20. Bp Buoyancy-driven flows, flow instabilities - 47.57. E- Suspensions, complex fluids - 47.55.P- 
Buoyancy-driven flows, convection 



■ 1 Introduction 

■ Rayleigh-Benard eonvection [l][2] is a prominent model 
' for many processes in geoscience, atmospheric dynamics, 
. and convection in technical applications [21I31[S] . Moreover 
' it is a classical lab experiment for studying generic phe- 
, nomena in nonlinear dynamics and pattern formation [51 
'El- Typically the Oberbeck-Boussinesq (OB) approxima- 
I tion is used for modeling thermal convection [TlfSlOlITU] . 
. In this approximation constant material parameters in- 
' dependent of the thermodynamic variables are assumed, 
, except for a temperature-dependent density in front of the 
I gravitational force, which is the essential driving force of 
I convection. 

An extension beyond the Boussincsq approximation is 
' accounting for temperature-dependent viscosity which is 
, used, for instance, for modeling various phenomena in the 

■ Earth's mantle [IIl[Il[Il[Il[fel[Ill[ni[IHl[Il[in]- A spa- 
I tially varying viscosity causes additional nonlinear cou- 
• plings between the temperature and the velocity field and 
' breaks the up-down symmetry of the flow field with re- 
spect to the center of the fluid layer. Often this symmetry 
breaking favors hexagonal patterns near the onset of con- 
vection. A linear or sinusoidal temperature dependence of 
the viscosity causes a reduction of the threshold compared 
to the case of constant viscosity [TTIIT^ITB] . For an expo- 
nential temperature dependence of the viscosity one can 
either have an enhancement or a reduction of the onset 
of convection, depending on the viscosity contrast, which 
measures the ratio of the viscosity at the boundaries [T^ • 
Strongly varying material properties in the Earth's mantle 
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are a major motivation for using a temperature-dependent 
viscosity in models of thermal convection in single com- 
ponent fluids [T51ll9ll20l[^l^ . Materials in the interior of 
the Earth are however mixtures of several substances. The 
dependence of the viscosity on the concentration of one 
component of the mixture and its influence on the onset 
of convection is one possible effect to be investigated. 

An interesting variant of the classical Rayleigh-Benard 
system represents thermal convection in binary-fluid mix- 
tures [7j. Here the concentration field of one of the two 
constituents enters the basic equations as an additional 
dynamic quantity. Given the possibility of an oscillatory 
onset of convection, this enriches the bifurcation scenario 
considerably and accordingly, the system has attracted 
wide attention by several reasons in the recent decades, in- 
cluding Non-Oberbeck-Boussinesq effects pl [23l[24l[25ll 2 6 1 
127] . In binary- fluid mixtures temperature gradients cause, 
via the Soret effect, concentration gradients and the con- 
centration field couples via the buoyancy term into the 
Navier-Stokes equations for the velocity field. In binary- 
fiuid mixtures (e. g. alcohol in water) concentration dif- 
fusion is roughly a factor of hundred slower than thermal 
diffusion. For finite amplitudes of the velocity field, this 
leads to concentration variations in a narrow boundary 
layer. 

In recent experiments Rayleigh-Benard convection was 
investigated in colloidal suspensions ^8,29,30] with a spe- 
cial focus on Soret driven convection (28tl29j and bistable 
heat transfer, involving new effects like e. g. sedimen- 
tation [5^. In theoretical studies of convection in col- 
loidal suspensions the model for binary-fluid convection 
is commonly used considering the colloidal nanoparticles 
as the second fluid constituent |31| . Since nanoparticles 
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are much larger than for instance alcohol molecules in wa- 
ter, their mass diffusion is more than two orders of mag- 
nitude smaller. Accordingly the Lewis number of suspen- 
sions reaches considerably smaller values, typically in the 
range of 10~* The Soret coefficient, representing the 
strength of the Soret effect, grows linearly with the parti- 
cle's size [32II33] and can therefore be changed in a wide 
range by varying mass density and size. 

It is well known that the viscosity of a suspension rises 
with increasing concentration of the solute particles [8l[34l 
[55] . The concentration of colloidal particles is sensitive to 
temperature gradients via the Soret effect. Hence varia- 
tions of the temperature cause spatial variations of the 
particle concentration. This leads to a spatially varying 
viscosity which is a function of the local concentration of 
the colloidal particles. The temperature dependence of the 
solvent's viscosity is neglected. The goal of this work is the 
analysis of the onset of convection in a colloidal suspen- 
sion in terms of a model for binary-fluid mixtures with a 
special emphasis on the effects of a particle-concentration- 
dependent viscosity. In Sec. [2] we briefly introduce the 
underlying equations of motion and the relation between 
viscosity and particle concentration. The linear stability 
analysis of the linear heat conducting state and the phase 
diagrams for the onset of convection are presented in in 
Sec. [21 Our results are discussed and summarized in Sec.[31 



2 Model 

Thermal convection in a horizontal layer of colloidal sus- 
pensions is described by the common mean field approach 
for binary-fluid mixtures 7 > 24. 25, 26, 27, 36 with one ex- 
tension: We take into account a linear dependence of the 
viscosity on the local concentration of colloidal particles. 

Assuming that the mass density of the colloidal par- 
ticles pc is similar to the mass density of the solvent ps, 
i. e. 



Pc 
Ps 



1 



(1) 



the difference between ps and the mean density of the 
suspension po is small and therefore sedimentation effects 
are ncghgible. In our case a spatially varying mass fraction 
of the coUoidal particles N{r, t) causes a spatially varying 
viscosity similar to Einstein's law [8]: 



= '/o ( 1 + e- [N - No) 



(2) 



Effects of higher order in N{r,t) are neglected [35l[57] . 
Here A^o represents the mean mass fraction of the col- 
loidal particles and 770 describes the mean viscosity of the 
suspension. 

The common set of basic transport equations for in- 
compressible binary-fluid mixtures, cf. Refs. [7)24,25,26^, 
[27], involves the density of the mixture p(r, t), the fluid ve- 
locity v(r, t), the temperature field T{r,t), the mass frac- 
tion of the particles A^(r, t), and the pressure field p{r, t): 



V-v = 0, 
(a* + V V) T = xV^T 

idt+vV)N = D 



\7^N + —V^T 
To 



(9t + v- V)v 



- — Vp + V • (i/V)v+ — g. 
Po Po 



(3a) 
(3b) 

(3c) 
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Eq. pap describes the incompressibility of the fluid, x in 
the heat equation (j3bp denotes the thermal diffusivity of 
the mixture and D in the diffusion equation (j3c|) the dif- 
fusion constant. Due to the size of the colloidal particles 
D takes much smaller values than in molecular fluid mix- 
tures. The dimensionless thermal-diffusion ratio fcy repre- 
senting the cross coupling between the temperature gradi- 
ent and the particle flux is related to the Soret coefficient 
St via fcj- = A^(l — N)TSt and can be either positive 
or negative. In the following fc^ — A''o(l — No)ToSt is re- 
garded as constant. In our model the kinematic viscosity 
v = rj/po in the Navier-Stokes-equations (|3d)) is a func- 
tion of space, as introduced by Eq. The gravity field 
g ~ —gSz is chosen parallel to the z-direction. 

For the local density p of the suspension we use a Hn- 
earized equation of state [TinU]. 



p = Po[l- a{T - To) + fi{N - A^o)] 



(4) 



with the thermal expansion coefficient a = —{1/ po)dp/dT 
and (3 = {l/po)dp/dN refiecting the density contrast be- 
tween the solvent and the suspended particles. According 
to the Boussinesq approximation this dependency of the 
density is taken into account only within the buoyancy 
term. The sign of (3 indicates whether the colloidal parti- 
cles have a higher or a lower mass density than the solvent. 
Here we assume /3 > 0, corresponding to e > 1. 

Boundary conditions. The fluid layer is conflned be- 
tween two impermeable and parallel plates at a distance d, 
and extends infinitely in the (x, y)-plane. The lower plate 
(z = —d/2) is kept at a higher temperature To + ST/2, the 
upper plate (2 = +d/2) at a lower temperature To~ST/2. 

Considering realistic no-slip boundary conditions for 
the flow field the boundary conditions at z = ±d/2 are: 



= = Vy = Vz = 9zWz , 

= d,N + ^ d,T, 
-to 



T = To T 



ST 



(5a) 
(5b) 

(5c) 



For geophysical applications free-slip boundary conditions 

= dzv^ = d^Vy = f 2 = dlvz , (6) 

for the flow fleld are often considered to be more realistic 
PTj . In the case of constant viscosity and free-slip bound- 
ary conditions an analytical determination of the onset of 
convection is possible [25]. This advantage is lost when 
introducing a concentration dependent viscosity and one 
has to rely on numerical methods. 
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Basic state. In the absence of convection (v = 0) the 
time- independent heat-conducting state is described by: 



To-ST- 



Ncond -^Nq-5N- with 5N 
a 



Tn 



ST. 



(7a) 
(7b) 



For further analysis it is convenient to separate the basic 
heat-conducting state in Eq. (jT]) from convective contribu- 
tions of the temperature and concentration fields as fol- 
lows: T{r,t) = T^ond{z)+Ti{r,t) andiV(r,0 = N^ond{z) + 
iVi(r,t). 

Since we have rotational symmetry in the fluid layer 
we can restrict our further analysis concerning the onset 
of convection to two spatial dimensions, namely to the 
(x, z)-plane. With this simplification the fluid velocity can 
be expressed by a stream function <P{x, z, t): 



(8) 



Subsequently all lengths are scaled by the distance d and 
time by the vertical thermal diffusion time cP /x- Scaling 
the temperature field T by {xvo)/{'^9d^)j the concentra- 
tion field N by —{hTXi^o) / (Taagcfi) and the stream func- 
tion by xd we are left with five dimensionless parame- 
ters: The Rayleigh number R, the Prandtl number P, the 
Lewis number L, and the separation ratio tpi 



P 



X ' 



L 



D 

1 

X 



R 



a 



gd^ 



Xvo 



ST, V 
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are well known from molecular binary-fluid mixtures |26l 
[27] , The fifth dimensionless quantity 



Q 



5 

■ e - 



Vox 



2 5/3^3 



(10) 



is introduced to characterize the spatially varying contri- 
bution to the viscosity of the suspension. For small par- 
ticles in water a reasonable value is Q ~ 0.01 and for 
glycerin Q — 10. A further illustration of its physical 
meaning is obtained by considering the viscosity contrast 
f] = rj{z = -\-^)lrj{z — — i) between the concentration- 
dependent viscosity at the upper and the lower boundary, 



l + \Ri}Q 
krST " 1 - ii??AQ ' 



5£ 

4To 



(11) 



which is essentially a function of the product of the three 
dimensionless control parameters: RipQ. In the following 
we will discuss our results essentially in dependence of "0 
and Q whereas P and L are regarded as constants. 

For further analysis we introduce a rescaled temper- 
ature deviation 9 = {R/ST)Ti and a rescaled concentra- 
tion deviation A^i = —{ToR/kTST)Ni in terms of these 
dimensionless quantities. In addition it is advantageous to 
introduce the combined function c ^ Ni — 9 instead of Ni . 



Suppressing for reasons of simplicity all the tildes we 
obtain: 

{dt - A)e ~ Rd,<P = {d,m - d^<Pd,)e , (12a) 
{dt - LA)c+Ae = id,m - d,^d,)c, (12b) 
{dt - PA)A$ - Pi^d^c - P{1 + i})dj 
-i^PQR [zA + 29^] A^ = {d,^d^ - d^$d,)A$ 

-^PQ [dl{c + e)dl<P + 2dA{c + e)d^d,'P 

+ dl{c + e)dl<i>]. (12c) 

No-slip, impermeable boundary conditions for the fields 
0, c, and ^ demand 



6» = a,c = = = at z ± 



1 



(13) 



while free-slip, permeable boundary conditions 



9 1 



(14) 



All results presented in this work base upon calculations 
with no-slip, impermeable boundary conditions. However, 
free-slip permeable boundary conditions lead to qualita- 
tively comparable behavior. 



3 Linear stability of the heat conducting 
state and the onset of convection. 

The critical values of the control parameters at the onset 
of convection are determined by investigating the growth 
dynamics of small perturbations 0{x,z,t), c{x,z,t), and 
(j){x,z,t) with respect to the heat conducting state given 
by Eqs. ([7]). The dynamics of such small perturbations 
is governed by the linear parts of the coupled nonlinear 
equations (HH. They are solved by the following expo- 
nential ansatz with respect to time t and the horizontal 
x-direction: 



{x,z,t) 
c{x,z,t) 
{x,z,t) 



with 



uo(z) 



Uo(z) e'«^e'^^ 

e{z) 
c{z) 
Mz)/{iq)^ 



(15) 
(16) 



With this ansatz and the eigenvalue a the linear parts 
of Eqs. (|12p are transformed into a boundary eigenvalue- 
problem with respect to z. The remaining linear ordinary 
differential equations (ODEs) are solved in this work by 
two different methods. 

The first one is a standard shooting method as de- 
scribed in detail for binary-fluid convection in Ref. |27j : 
The resulting coupled ODEs for Uo(2;) are integrated by 
a standard solver for ODEs. For every set of initial condi- 
tions a boundary determinant /(ct, i?, g, Q, P, L, ip) is cal- 
culated. During an iteration procedure R or a are varied 
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Fig. 1. Neutral curves Ro{q) obtained according to Eq. H17|l 
are shown for different values of the viscosity parameter, Q = 
(solid line), Q = 5 (dashed line) and Q = 10 (dotted line). In 
part a) for the positive separation ratio ip — and in b) for 
-iP = 10"^ 

such that / vanishes. This yields a and i? as a function of 
the remaining parameters. 

The second possibility is the so-called Galerkin method 
where the components of uo{z) are expanded with re- 
spect to a suitable chosen set of functions which fulfill the 
boundary conditions, either Eq. or Eq. ([Tl)) - see, e. g., 
Ref. [5^liniliT] . Here the linear ODEs are transformed into 
a generalized algebraic eigenvalue problem which can be 
solved numerically. 

At the onset of convection small perturbations 9, c, 
and (f) with respect to the Hnear ground state neither grow 
nor decay, i. e. we are interested in the neutral stability 
condition for the fastest growing mode (eigenvalue a with 
the largest real part), 

Re(cr) = with (J = a{R,q,Q,P,L,iP), (17) 

for the determination of the threshold. Depending on pa- 
rameters the system may show a stationary bifurcation 
with a vanishing imaginary part of the eigenvalue, lui{a) = 
w = 0, or a Hopf bifurcation with a finite Hopf frequency 
Im(cr) = ±u}^Q. 

For a given set of material parameters the neutral sta- 
bility condition (jl7p provides the Rayleigh number Ra{q) 



as a function of q, the so-called neutral curve 

Ro{q) = Ro{q,Q.P,L,^). (18) 

Similarly we obtain from Eq. (jl7l) luq = u)f)[q^Q,P,L,ip) 
in the case of a Hopf bifurcation. 

In this work wc choose a Prandtl number of P = 10. 
Since the mass diffusion of the colloidal particles is much 
lower than the thermal diffusion, a Lewis number of L = 
10~^ is a reasonable assumption. For small particles in 
water Q ~ 0.01 and in Glycerin Q — 10 are appropriate 
values. As the main changes for finite values of Q occur 
if ij] is comparable to L we concentrate on rather small 
values of ■0. 

Neutral curves for the present system are given in 
Fig. [1] for two different positive values of the separation 
ratio -0, namely for %p = 10~* and -0 = 10~^, and for three 
different values of the viscosity parameter Q = 0, 5, 10. 
The onset of convection is given by the minimum of the 
neutral curves Rc = Ro{qc) = mini?o('z) which can be 
compared with its corresponding value in the case of a 
simple Newtonian fluid, R^^ ~ 1708, and molecular bi- 
nary mixtures |27j . 

The reduction of the threshold in the range of > is 
well known from molecular binary fluids |25ll27j and the 
reason is as follows. Due to thermal expansion a fluid layer 
heated from below becomes lighter at the lower boundary 
and there is a temperature induced density gradient. Ac- 
cording to a positive Soret effect (0 > 0) the colloidal 
particles tend to move to the upper plate which is colder. 
As the colloids are heavier than the molecules of the sol- 
vent the Soret induced concentration gradient amplifies 
the temperature induced density gradient. Hence the fluid 
layer becomes unstable already at temperature differences 
much lower than for the onset of convection in a simple 
fluid. 

For positive values of "0, shown in Fig. [1] the bifurca- 
tion from the heat conducting state is stationary. While 
the minimum of the neutral curve is located at a finite 
wavenumber qc > for values of ^ that are small com- 
pared to L (cf. Fig. [TJa)), the critical wavenumber moves 
to = for larger values of ip (cf. Fig.[IJb)). Finite values 
of Q lead to a lowering of the neutral curves whereby the 
critical wavenumber is nearly unaffected. The differences 
between the critical wavenumbers qc at the minimum of 
the neutral curves in Fig.[TJi) is less than 0.1%. 

Fig. [5] shows the threshold reduction as a function of 
Q for three different positive values of 0. Here we plot the 
critical Rayleigh number Rc{Q) normalized by the critical 
Rayleigh number Rc{Q = 0) which we call the relative 
critical Rayleigh number 

The reduction of Rc,rei{Q) as a function of Q shows a 
parabola-like behavior as shown in Fig. [21 which is ac- 
companied by a tiny variation of the critical wavenumber 
qc along each curve. 

For finite values of the ground state displays a viscos- 
ity contrast fj ^ as described by Eq. ([TT|) . Fig. |3] shows 
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Fig. 2. The relative critical Rayleigh number Rc,rei = 
Rc{Q) / Rc{Q ~ 0) is shown as a function of Q for three differ- 
ent values of the separation ratio ip: For ip = lO^^ < L (solid) 
one has Qc — 2.97, for ip = lO^'* — L (dashed) qc = 2.10 and 
for = 10"^ > L (dotted) Qc = 0.00. 
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Fig. 4. The relative critical Rayleigh number Rc,rei = 
Rc{Q) / Rc{Q ~ 0) is shown as function of the product ~tpQ for 
two different values of the separation ratio ijj: For ijj = — IQ-" 
(solid) and for ip — —10-^ (dashed). The critical wavenumber 
along both curves takes the value qc — 3.116. 




Q " ^ ^ " Q 

(a) V = 10-" (b) iP = -10-" 

Fig. 3. In part (a) the viscosity contrast fj, given by Eq. ()11|) . 
is plotted for = 10"" and in part b) the inverse expression 
j;-" is plotted for ^ — -10"". In both cases the Rayleigh 
number R entering Eq. pip was taken at the respective onset 
of convection. 



fi{Q) as a function of Q for two different values of ip at the 
onset of convection R ~ R^. The Q-range in Fig. |3] was 
limited to the range \RQ'4>\ < 1/4 where fcy = const, is a 
reasonable assumption. 

Another interesting question is how the threshold re- 
duction behaves as function of the viscosity contrast fj. 
Within our model the viscosity depends linearly on the 
concentration, according to Eq. Since the concentra- 
tion, via Eq. (|7b|) . is a linear function of the temperature, 
the viscosity is effectively a linear function of the tem- 
perature too. For a given viscosity contrast the thresh- 
old reduction of stationary convection {if) > 0) is similar 
to that in a previous study of a single component fluid 
with a direct linear temperature dependence of the viscos- 
ity [12]. However, there are important differences between 
this single component study and the present model: Here 
the viscosity contrast is caused by a well known physical 
mechanism, the Soret effect. The viscosity contrast is ac- 
companied by heterogeneous distributions of the colloidal 
particles, which can be of relevance in the context of re- 
lated convection phenomena in geoscience. For a single 
component fluid in Ref. |16j one has only a stationary on- 
set of convection. In contrast the concentration field for 
the colloidal particles represents an additional degree of 



freedom which, in the range of V' < 0, similar as in molec- 
ular binary- fluid mixtures p51l27] . leads to an oscillatory 
onset of convection via a Hopf bifurcation. 

Therefore arises a further interesting question: How 
does a flnite viscosity contrast affect an oscillatory bifur- 
cation? We find in the range of '0 < 0, as in the range 
of > 0, a lowering of the neutral curve with increasing 
values of Q. The critical wavenumber as well as the 
Hopf frequency Wc, show only small and negligible changes 
when varying Q. The Q dependence of the relative criti- 
cal Rayleigh number Rc,rei{Q) is slightly different from its 
behavior in the positive range of ■0. For > the criti- 
cal Rayleigh number i?c('0) decays strongly as a function 
of ip, whereas the threshold increases only moderately for 
decreasing ip in the range ip < Q [251l27j . Therefore, the 
viscosity contrast at threshold given by Eq. ([TT]) is for 
< essentially a function of the product -0(5. 

The relative threshold Rc,rei{Qi') is shown in Fig.Uas 
a function of the product ~Qip for two different values of 
"0 < 0. The two curves are nearly identical as expected. 
Therefore, the changes of the threshold of the oscillatory 
convection are again essentially determined by the viscos- 
ity contrast fj given by Eq. (jlip in the restricted range 
\R'il}Q\ < 1/4. 

Another view on the parameter dependence of the thresh- 
old is given in Fig. [SJ where the relative critical Rayleigh 
number i?c,rei('0) is shown as a function of the separation 
ratio "0 for three different values of the viscosity param- 
eter Q = 3.0, 5.0, and 10.0. In the limit of a vanishing 
separation ratio ip the particle concentration and there- 
fore the viscosity contrast fj does not change by varying 
Q, so we have Rc,rei — 1- For increasing values of tp the 
threshold drops down and finally reaches a constant value, 
depending on Q. The strongest decay can be observed in 
the region ip L accompanied by the critical wavenumber 
qc going down, qc finally vanishes at nearly the same value 
of ~ 4 X 10-" independently of Q as shown by Fig. [SJd). 

For '0 > an increasing viscosity contrast leads to a 
decrease of the velocity field at the upper boundary as the 
fiuid motion tends to concentrate in the region of lower 
viscosity at the lower boundary. Therefore the center of 
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Fig. 5. Part a) shows the relative critical Rayleigh number 
Rc,rei{'4'jQ) = Rci'4' J Q) / Rci4' J Q = O) as a function of the 
separation ratio il) for three different values of Q: Q = 3.0 
(sohd), Q = 5.0 (dashed), and Q = 10.0 (dotted). Part b) 
shows the corresponding critical wavenumber qc{'4')- 



the convection rolls at zq is slightly shifted below the cen- 
ter of the convection cell at z = as shown in Fig. [5^) 
for -0 = 10~^. This shift zq is more clearly indicated in 
Fig. [6)3), which shows the horizontal component of the 
convection velocity, Vx[z), as function of z. The shift Zq 
to negative values increases with increasing values of Q as 
shown in Fig. [51:). This behavior is common for all positive 
values of ^ with a finite value of the critical wavenumber 
qc- For t/i < 0, in the range of an oscillatory Hopf bifur- 
cation, the viscosity is higher at the lower boundary and 
therefore the magnitude of the velocity is reduced near the 
lower plate. Accordingly, the center of the convection rolls, 
zq, is elevated with increasing values of ijjQ to zq > 0. Our 
analysis shows that the shift |zo| depends essentially on 
the viscosity contrast 77. 

Assuming free-slip boundary conditions p4p instead of 
no-slip boundary conditions p3p . wc obtain very similar 
results. For finite values of Q the threshold is reduced 
compared to its value at Q = and the center of the rolls 
is shifted to the region of lower viscosity. 



4 Conclusions 

The onset of thermal convection in a colloidal suspen- 
sion was investigated in terms of a generalized continuum 
model for binary- fluid mixtures going beyond the common 
Boussincsq-approximation by taking into account a linear 
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Fig. 6. The flow velocity field at the onset of stationary con- 
vection is shown in part a) for the parameters -0 — 10~* and 
Q — 10. Part b) depicts the variation of the velocity component 
Vx (z) for the same parameters as a function of the vertical co- 
ordinate z and part c) shows the shift of the position zq where 
the horizontal velocity vanishes, i. e. Vx{zo) = 0, as a function 
of Q. 



dependency of the viscosity on the local concentration of 
the colloidal particles. 

Spatial changes of the concentration of colloidal par- 
ticles in the convection cell arc induced by temperature 
variations via the Soret effect similar to that in molecular 
binary-fluid mixtures. The influence of the particle con- 
centration on the viscosity is described by a dimensionless 
viscosity parameter Q, introduced in this work. For finite 
values of Q the threshold is reduced compared to the limit 
(5 = 0, for both the stationary bifurcation ioi ip > and 
the oscillatory Hopf bifurcation for '0 < 0. This reduction 
depends essentially on the induced viscosity contrast fj at 
the choosen values of Q and ip. We would like to empha- 
size that \RipQ\ may not be too large so that the linear 
approximation of the concentration variation in Eq. (7b) 
remains a reasonable assumption. 
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For both bifurcations the center of the convection rolls 
is shifted slightly away from the center of the convection 
cell: For ip > the convection rolls are shifted below and 
ior ip < beyond the center of the cell, as the fluid motion 
always tends to concentrate in the region of lower viscos- 
ity. At the onset of convection the changes of the critical 
wavenumbcr of the convection rolls as a function of the 
dimensionless parameter Q are rather small. 

In a previous study on convection for a single compo- 
nent fluid in Ref. [16] a linear dependence of the viscosity 
on the temperature was assumed and its influence on the 
onset of stationary convection was investigated. Assuming 
the same viscosity contrast a similar relative reduction in 
the threshold is obtained as in our model. However, our 
model does not only describe a spatial variation of the ma- 
terial parameters. It describes also a mechanism leading, 
via a heterogeneous concentration distribution, to a spa- 
tial dependence of the viscosity. This interrelation between 
concentration and viscosity variations may be of relevance 
to related phenomena in geophysical contexts. Addition- 
ally, the further degree of freedom in our model, intro- 
duced by the concentration dynamics, provides a richer 
bifurcation scenario and we can investigate the effects of 
a spatially varying viscosity on a stationary as well as on 
a Hopf bifurcation. 

In this work the viscosity fjQ of the solvent of the col- 
loidal particles was assumed to be constant. A natural 
extension of the present model would be to take also into 
account the effects of a temperature-dependent viscosity 
of the carrier fluid. In the range of larger values oi ip > 
with small temperature differences at the threshold, this 
effect may be small. However, in the range of small and 
negative values of -0 this effect may play a role and one 
may then have a spatially varying mobility of the colloidal 
particles corresponding to a temperature and spatially de- 
pendent Lewis number. A temperature dependent Lewis 
number is also expected for thermosensitive microgel sus- 
pensions [42ll43] in which core-shell colloids change their 
size as a function of the temperature and their mobility 
as well as the viscosity of the suspension. This combined 
effect on the viscosity as well as the effects of nonlinear 
variations of the particle concentration will be discussed 
elsewhere. 

The vertical viscosity variation, caused by a temper- 
ature gradient and leading to spatially dependent viscos- 
ity, breaks the up-down symmetry in the convection cell. 
Therefore the interesting question arises: In which param- 
eter range are hexagonal patterns favored for this binary- 
fluid model in the weakly nonlinear regime? 

Stimulating discussions with F. H. Busse, M. Evonuk, G. Fre- 
und, W. Kohler, W. Pesch, I. Rehberg, and W. Schopf are 
appreciated. This work has been supported by the German 
science foundation through the research unit FOR 608. 
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